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I Abstract. Strong constraints on the cosmic star formation history (SFH) 

have recently been estabhshed using ultraviolet and far- infrared measurements, 
refining the results of numerous measurements over the past decade. Taken 
^ , together, the most recent and robust data indicate a compellingly consistent 

O ■ picture of the SFH out to redshift z « 6, with especially tight constraints for z < 

' 1. There have also been a number of dedicated efforts to measure or constrain 

the SFH at z « 6 and beyond. It is also possible to constrain the normalisation 
of the SFH using a combination of electron antineutrino flux limits from Supcr- 
Kamiokande measurements and supernova rate density measurements. This 
review presents the latest compilation of SFH measurements, and summarises 
' the corresponding evolution for stellar and metal mass densities, and supernova 

QQ I rate densities. The constraints on the normalisation of the cosmic SFH, arising 

, from the combination of the supernova rate measurements and the measurement 

limit on the supernova electron antineutrino flux, are also discussed. 

o 

:^ 

Qh' !• Introduction 

O 

I In the past few years measurement of the evolution of galaxy luminosity func- 

^ . tions at a broad range of wavelengths has rapidly matured. One consequence 

of this has been the refinement in our understa nding of how the space density 
of tli e galaxy star formation rate (SFR) evolves (iLillv et al.lll996l : iMadau et al.l 



> 



X 



19961 ). In particular the cosmic star formation history (SFH) is now quite tightly 



constrained (to within « 30 — 50%) up to redshifts of z ~ 1. Combined with 
\ measurements at higher redshifts th e SFH is well d etermined (within a factor of 

about 3 at z ^ 1) up to z ~ 6 fe.g.. lHoDki^l2004l l. 



Additional results from the Super-Kamiokande (SK) particle detector pro- 
vide a strong limit on the electron antineutrino {Ve) flux, 1.2cm~^s~^ (for 
El, > 19 .3 Me V), originati ng from supernova type II events associated with 
the SFH (Mal ek et al.ll20()3 ). This limit on the diffuse su pernova neutrino back- 



ground (DSNB) con strains the normalisation of the SFH (jFukugita &: Kawasaki 



20031 : IStrigari et al.l [2005). An exploration of quantities predicted from the 



SFH, the stellar and metal mass density evolution, and supernova (SN) rate 
evolution, provides further insight into the allowable normalisation of the SFH 
( Strigari et al.ll2005l 1. This series of interconnected physical properties of galax- 



ies and SNe provides an emerging opportunity for determining the level of the 
SFH normalisation, and the SFH measurements particularly for z ^ 1 now have 
the precision to allow this exploration of their accuracy. Constraining the nor- 
malisation of the SFH will support a range of quantitati ve analyses of galax; 



evolution, including the mass-dependence of the SFH (e.g.. lPapovich et al.ll200i 



1; 
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Figure 1. Evolution of SFR density with redshift (scaled assuming the SalA 
IMF). Circles are from the compila tion of Hopkins (2004). The hatched region 
is the 24 iim SFH fromlLe Floe h et all (|2005l) . Triangles are 24 ^m data from 
iPerez-Gonzalez et all (I2005D . The open star at z = 0.05 is based on 1.4 GHz 
data from Mauchl ( 20051) . The filled circle at z = 0.0 1 is the Ha e s timat e 
from'Hanish et al.' ("20061). Squares arc UV data from 'Baldrv et al.' (2005); 
Wolfetal. (2003); Arnout s et al. (2005); Bouwens et al. (2003a,b, 2005a); 
Bunker et al.l (120041): lOuchi et all p004i ). Crosses are the UDF estimates 
from [Thompson et al.l (2006[ ) 



Juneau et al. 2005: Heavens et al. 


2004 ) , and the reasons 


in the SFH to low redshifts (e.g., 


Bell et al. 20051). 



The analysis o f the c onstraints on the SFH normalisation are detailed in 
iHopk ins &: Beacoml (|2006l ). and here the main arguments are summarised. Since 
optical SNH can be hidden from observations by dust obscuration, the present 
SNH rate density measurements may merely be lower limits, and serve as a 
lower bound on the allowable SFH normalisation. In contrast, since neutrinos 
are unaffected by dust, the DSNB provides an absolute upper limit on the true 
SNH rate. We assume Hq = 70kms~i Mpc'^, VLm = 0.3, f^A = 0.7. 



2. The Data 



The compilation of HopkinsI ( 20041 ) was taken as the starting point for this analy- 
sis, and uses their "common" obscuration correction where n ecessary. Additional 
measu rements are indicated in Figure [U and are detailed in iHopkins &: Beacoml 
(l2006l ). 



2.1. Dust Obscuration Corrections 

To implement effective obscuration corrections for the UV measurements at 
z ^ 1, we take advaii t age o f the well-established FIR SFR densities up to z = 1 



from lLe Floc'h et al.l (|2005l ). The UV data at z < 1 are "obscuration corrected" 



by adding the FIR SFR density from lLe Floc'h etlil] (l2005l ) to each point. As 
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Figure 2. SFH data used in the parametric fitting, (a) Assumes SalA IMF. 
(b) Assumes BG IMF. The shape of the fits is determined from the SFH 
data alone, and a scaling factor is fit to ensure consistency with the SK 
limit Solid lines assume a Fg temperature of T = 4 MeV or T = 6 MeV, and 
dotted r = 8 MeV. The grey shaded and hatched regions are the 1 a and 3 cr 
confidence regions around the T = 4 MeV fits respectively. 



shown by iBelll ( 2003al ) for individual systems, this technique results in p^, esti- 
mates consistent with the obscuration c orrected Ha me a surem ents. This result 
is consistent with the interpretation of Takeuchi et al.l ( 20051 ) that about half 
the SFR density in the local universe is obscured by dust, increasing to about 
80% by z 1. For obscuration corrections to the UV data between 1 < z < 3 
we rely on the fact that the FIR measurements of iPerez- Gonzalez etliLl (|2005l ) 
are quite fiat in this domain, as well as being highly consistent with those of 
Le Floc'h et ah J 20051) at z < 1, and add the constant SFR density correspond- 
ing to that of Le Floc'h et al. ( 20051 ) at z = 1. This is also consistent with the 
recent me a surein ents of obscuration corrections for UV luminosities at z ~ 2 by 
lErb et all ( 20061 ) . who find a typical correction factor of w 4.5. At higher red- 



sh ifts we apply a "common" obscuration correction to the UV data as detailed 
in Hopkins (2004). 



2.2. IMF Assumptions 

While uncertainties in SFR calibration act to increase the scatter in the SFH, 
and uncertainties in dust obscuration can raise it to greater or lesser degrees, the 
choice of IMF is the only a ssumption that can systematically decrease the SFH 
normalisation. A modified ISalpetei ( 19551) IMF, with a turno ver below IMq, 



remains a reasonable model (e.g., Baldrv &: Glazebrookl [2003 ) , and other cur- 



rently favoured IMFs include those of Kroupa ( 200ll ) and Baldrv &: Glazebrook 



(2003). The factor to c onvert SFH measurern e nts fr om a traditional ISalpetei 
1955. ) IMF to the IMF of lBaldrv Glazebrookl (|2003l . her eafter BG IMF) is 0.50 



(—0.3 05 dex). To convert to the modified Salpeter A IMF ( Baldrv &: Glazebrookl 
20031 . hereafter SalA IMF) is a factor of 0.77 (-0.114dex). The iKroupai (l200llT 
IMF and the modified Salpeter B IMF tealdrv fc Glazebrookl |2003| . hereafter 
SalB IMF) have scale factors intermediate between these choices. With these 
two extreme IMF choices we expect to provide bounds encompassing the result 
from choosing any reasonable IMF in our subsequent analysis. 
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Figure 3. SFH measurements emphasising recent ^ ^ 6 estimates (assum- 
ing the SalA IMF). The filled circles and grey shaded and hatched regions 
are as in the previous Figure (note the different axes scales). The dashed 
l ine corres ponds to the level of (unobscured) SFH required f or reionisation 
()Madau et al. 199 ^. Triangles: Lyg emit ters ()Ive et al.l 120061) : Stars: LBGs 
i n GOODS (iBouwe ns fc Illingworthl l2n"(i6h: Filled squares: LB Gs in GOODS 
(iMannucci et all 12006): Diamond: LBGs in the HUDF ( Bouwens et al.l 
Hexagon: LB Gs (A. Verma et al., 2006, this proceedings); Open 
squares: lensed LBGs ( Richard et al . 2006). In all cases where the same sym- 
bols appear at the same redshift, the upper values correspond to the maximum 
expected SFH after integration over the full luminosity function. 



3. SFH Fitting 

To derive a Ve flux from the DSNB for comparison with the limits fr om SK, we 



first fi t a functional form to the SFH. We use the parametric form of I Cole et al 



(|200ll ): p^ = {a + bz)h/{l + {z/cf ), here with h = 0.7. The individual p* mea- 



surements chosen to constrain this fit are also important since the resulting fit 
will obvious ly vary depending on the d ata used. The details of the data selected 
are given in Hopkins &: BeacomI ( 20061 ) . and the results shown in Figure [2l The 



final parametric fitting is a simple fit to the 58 selected /j* measurements 
spanning < z ^ 6. 

3.1. UV Data at High Redshift 

As a brief aside, it is worth discussing star formation at z ^ 6, as there are now a 
large number of estimates in the literature spanning this range. These primarily 
use the photometric dropout technique to select Lyman break galaxies (LBGs), 
often using relatively small, very deep fields (including the Great Observatories 
Origins Deep Survey, GOODS, and the Hubble Ultra Deep Field, HUDF). Mea- 
surements of Lya emission typically show SFH measurements significantly lower 
than those from LBGs. It is not clear whether the same populations are being 
probed here, or the full extent of the selection biases that may be present. 
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Figure 4. The Ve energy spectra predicted for the various SFH fits and 
temperature assumptions. The solid and dashed curves are the SalA IMF and 
BG IMF respectively. The T = 4 MeV and T = 6 MeV curves are consistent 
with the SK Ve limit. The T — 8 MeV curves are inconsistent with the Ve 
limit, and indicate the shape of the Ve spectrum derived by assuming the 
parametric form for the SFH corresponding to our best fit (T = 4 MeV) , and 
setting the Fe temperature to the higher value. The thin vertical line marks 
E = 19.3 MeV, above which the X/g contribute to the SK limit. 



The ^ 6 SFH contributes a negligible fraction to the integrated SFH that 
produces the DSNB, and is neglected in the subsequent analysis. It is, however, 
of great significance to the epoch of reionisation, and star formation may in fact 
play the dominant role here. Figure O shows the measurements of the z ^ 6 SFH, 
based on rest-frame UV luminosities not corrected for obscuration effects. The 
dashed line in this Figure is the le vel of SFR density required for reionisation at a 
given redshift (|Madau et al.lll999l 1 assuming /esc = 1 , appropriate for unobscured 



UV emission, and a clumping factor C = 30. 

Although the issue of dust obscuration at high redshift is still highly un- 
certain, some data is begin n ing to be obtained. In addition to the E{B — V) 
estimates from Ouchi et al. ( 20041 ). intriguing evidence for significant obscura- 



tion (Ay ~ 1 mag) at z = 6.56 has recently been establis hed through Spitze r 
observations of a lensed Lyman a (Lya) emitting source (jCharv et aDl2005l ). 



This implies that the first epoch of star formation in this source must have oc- 
curred around z 20. This is a lso supported by spectroscopic Lya emission 
measurements of LBGs at z ^ 5 (jAndo et al.ll2005l ). suggesting that the bright 



LBGs (at least) lie in dusty, chemically evolved systems at this redshift. One 
promising suggestion here is that the cosmic SFH beyond z ^ 6 can be probed 
through galaxy archeology, i.e., by determining the star formation histories of 
the z ^ 6 galaxy population (see also R.-R. Chary 2006, this proceedings). 

3.2. The Diffuse Supernova Neutrino Background 

The DSNB derived from the cosmic SFH is calculated as follows. The p*{z) 
is first converted to a type II supernova rate history, /0sNii(-2)- For the IMFs 
explored here psmi{z) = (O.OI32/M0) /j*(z) for the BG IMF and psNii(^;) = 
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Figure 5. Left: Evolution of stellar mass density buildup inferred from the 
SFH. Ri ght: Evolution of metal m ass density buildup inferred from the SFH. 
Refer to iHopkins fc BeacomI ()2006l ) for references to the measurements shown 
in both panels. The grey shaded and hatched regions come from the 1 a and 
3 a confidence regions around the T = 4 MeV fits to the SFH respectively. 
The Sal A IMF is assumed in all cases. 



{0 .00915 / Mq) p^{z) for the SalA IMF. The predicted differential neutrino flux 

(per unit energy) is then calculated by integrating psmiiz) multiplied by the Fp 

emiss ion per supernova, appropriately redshifted, over cosmic time (see lHopkins &i Beacom 
hood , for details). Finally, this energy spectrum is integrated above 19.3 MeV to 



estab lish the Ue flux for comparison with the SK limit of 1.2 cm ^ s ^ (jMaIek et al 



20031 ). We explore the implication of assuming a temperature of T ~ 4 MeV, 
6 MeV, or 8 MeV. Given the flux for each temperature, we simply scale the 
best fitting SFH to ensure the SK limit is not violated. The best fit SFHs inde- 
pendent of the I7g limit are identical to t he fits constra i ned b y a Ve temperature 
of T = 4 MeV. As found and discussed in lYiiksel et al.l ( 20051 ). our results favour 
effective temperatures at the lower end o f the predicted range. 

In addition to the ICole et al.l (|200ll ') parameterisation, we also explored a 
piecewise linear SF H model in l og(l + z) — log(/9*) space, in order to test the 
possibility that the ICole et al.l (|200ll ) parametric model could be biasing the 
shape of the resulting SFH fit in some way. In this model we allow the following 
six parameters to vary: The z = intercept, the slopes of three linear segments 
and the two redshift values at which the slopes change. 



4. Results and Discussion 

Figure [J shows the current SFH data compilation (assuming the SalA IMF) 
emphasi sing the additio nal data used in t his analysis comp ared to the compi- 
lation of lHopkinsI (|2004l '). The best fitting ICole et al.l (|200lh form for this IMF 



is also shown as is the best-fitting piecewise linear fit. Figure [2] shows the data 
used in the fitting and the best fits assuming three temperature values for the 
Vp population for e ach IMF assumed. The fitting parameters are tabulated in 
Hopkins fc BeacomI §006). 



For both assumed IMFs it can clearly be seen that the assumption of T = 
8 MeV, when the SFH is required to be consistent with the Fg flux limit, is 
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Figure 6. Evolution of SN rates inferred from the SFH. The upper curves 
correspond to the prediction s for psNiii and the l ower for psma, assuming a 
delay time tia = 3 Gyr. See iHopkins fc BeacomI (|2006|) for references to the 
measurements shown. 



inconsistent with the SFH measurements. Also for both IMFs, the best fitting 
SFH assuming T = 6MeV is identical to that assuming T = 4MeV. This can 
be understood by considering the SalA IMF, for example, with the higher SFH 
normalisation, but which also has a lower conversion factor between and /osNii, 
causing the predicted flux to be within the SK limit even with the assumption 
of the slightly higher neutrino temperature. For both IMF assumptions we 
determine the 1 a (grey-shaded) and 3 a (hatched) confidence regions around 
the best fitting SFH (corresponding to T = 4 or 6MeV). Subsequent Figures 
reproduce these confidence regions in the predictions for stellar and metal mass 
density evolution {p*{z) and pziz), respectively) and SN rate evolution (psN (-?))■ 
The shape of the Vg energy spectra corresponding to the various fits for the 
different IMF and temp erature assumptions are shown in Figure H] (see also 



Beacom Strigarill2006l ) 



4.1. Stellar and Metal Mass Densities 



Figure [5] shows the evolution of the stellar mass de nsity, p*(z), (a far more exten- 
sive compilation of data appears in Figure 4 of iFardal et al. 20061) . and the metal 



mass density, pz{z), as inferred from the SFH (jPei &: Falll Il995l : iMadau et al 



19961 ). The predictions from the best fitting SFH are also shown. For the evolu 



tion of p*{z) we need to know the fraction of the stellar mass recycled into the 
interstellar medium as stellar winds o r SN ejecta, R, corresponding to each IMF 
(|Cole et al.ll2"00ll : iMadau et allll998l ;i. We find R = 0.40 for the SalA IMF, and 



R = 0.56 for the BG IMF. Th e stellar mass inf erred is then a fraction (1 — R) of 
the time integral of the SFH (jCole et al.ll200ll ). Converting the observed stellar 
mass density measurements (where a Salpeter IMF is assumed) to our assumed 
IMFs is achieved by scaling by the product of the SFR conversion factor and 
the ratio of the 1 — R factor for the chosen IMF to that of the Salpeter IMF 
(where 1 — R = 0.72). The p*{z) measurements in Figure [5] clearly lie sys- 
tematically below the predictions from the SFH. Simulations suggest that the 
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measurements might be und erestimating th e total p^(z) (e.g . . iNagamine et al. 
20041 : ISomervihe et alJ l200lh . Discussion by iDickinson etalT (|2003l ') imply that 
it is perhaps not unreasonable to expect about a factor of two larger stellar 
mass densities as a result of reasonable obscuration levels. Other issues that 
have also been raised include incomplete galaxy population sampling and cos- 
m ic variance that may a ffect surveys probing small fields of view (see discussion 
in Nagamine et al. 20041 ) . At low redshift, the discrepancy between the measure- 
ments of p^(z) and the SFH predict ion is more of a concern, and is discussed 



ments or p^(z} ana tne b-bM predict ion is more oi a concern 
further bv lHopkins BeacomI ( 20061 . and references therein). 



To determin e pziz) from the SFH, we assume that p^. = 63.7 pz (e.g., 
Conti et al.l |2003| ) . At z = the compilation of data from ICalura &: Matteucci 
(120041 ') is shown, and these authors favour a valu e of 1.31 x 10^ MqMpc ^. Values 
at z = and z = 2.5 from Dunne et al.l ( 20031 ) are also shown, suggesting that 
the evolution in pz from the SFH may be consistent with that estimated from 
the dusty subniillini eter galaxy (SMG) population, although recent results from 
Bouche et all (|2005l ) indicate that the SMGs may contribute much less to the 
metal mass density at high redshift. Investigating the predictions of pz{z) from 
the SFH are complicated by the limited number of estimates for this quantity 
at z > 0. This is observationally a difficult measurement to make, particularly 
as much of the metals may exist in an ionised in tergalactic medium component 
(see, e.g.. lDunne et al.ll2003l : feouche et al.ll2005l l. 



4.2. Supernovae Type la and Type II 

Figure [6] shows the evolution in the SN rate for both types la and II SNe. 
The SNII rate density, psNiij can be calculated directly from the SFH. The 
SNIa rate density, psNia, involves more assumptions about the p roperties of 
SNIa events than in the case of SNII. See lHopkins BeacomI (|2006l ) for details 
of these calculations. The psNii measurements provide a strong lower bound 
on the normalisation of the SFH. Particularly given that uncertainty regarding 
obscuration corrections is more likely to raise than lower the psNii measurements, 
the SFH normalisation cannot realistically be much lower than that obtained 
from assuming the BG IMF (Figure [2j)) . Moreover, the psmi measurements are 
unlikely to be affected by sufficient obscuration to sup port an SFH normalisatio n 
much higher than that obtained with the SalA IMF (jHopkins fc Beacom 20061 ). 

The prediction for psm-a. from the SFH is also particularly intriguing. The 
assumption of the fixed tja = 3 Gyr has the effect of matching the z ^ 3 turnover 
in the fitted SFH with the apparent decline in psNTa seen i n the highest redshift 



measurement from the GOODS sample of Dahlen et al. ( 20041 ) . It is possible 
indeed probable, that this is simply a coincidence as it is a single psNia measure- 
ment, with large uncertainties, that is suggestive of the decline, an d the turnover 
i n the SFH is driven almost entirely by the z ~ 6 measurement of I Bunker et al 



( 20041 ) . It is thus still highly possibl e that the decline in both t he SFH and psNia 
lie at somewhat higher redshift (see lHopkins &: Beacom 20061 ). 



5. Summary 



The SFH compilation of Hopkins ( 20041 ) has been updated, emphasising the 
strong constraints from recent UV and FIR measurements, and refining the 
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results of numerous measurements over the past decade. The results of para- 
metric fits to the SFH, constrained by the SK I7g limit, suggest that the preferred 
IMF should pro duce normalisations within t he range of tho se from the modified 
Salpet er A IMF (iBaldrv fc Glazebrook|[2003 ) and the IMF of lBaldrv Glazebrookl 



( 20031 ) ■ They also suggest that lower temperatures (T = 4— 6 MeV) are preferred 



for the Fe population. 

Based on the fits to the SFH we predict the evolution of , pz and psN , and 
compare with observations. The psNii measurements provide a key constraint 
on the SFH normalisation, and correspondingly on the favoured IMF. In partic- 
ular, these data bound the SFH from below, while the DSNB bounds the SFH 
from above. Together, these provide a novel technique for testing or verifying 
measurements of a universal IMF. More measurements of psN for both type II 
and la SNe over a broader redshift range would help to more strongly constrain 
both the preferred universal IMF and the properties of SNe. Observing the high 
redshift turnover in the SNIa rate would also have strong implications for the 
location of the expected high redshift turnover in the SFH. 

With the best fitting SFH models explored here, the predictions for the 
DSNB appear to lie excitingly close to the measured Fg flux limit. Direct obser- 
vation of the DSNB will clearly allow much greater insight into the physics and 
astrophysics of star formation and supernovae. Already the DSNB constraint 
indicates a preferred IMF range and normalisation for the SFH. It also illustrates 
that stronger constraints on the SFH have implications for understanding the 
details of both SNII and SNIa production, and the physical basis of neutrino gen- 
eration by SNII is intimately associated with all these predictions. Being able 
to detect the DSNB and its energy spectrum will allow a more sophisticated 
analysis of the detailed connections between all these aspects of star formation 
and the cosmic SFH. Methods for increasing the sensitivity of particle detec- 
tors to DSNB antineutrin os and neutrinos have been detailed elsewhere, and 
Hopkins k BeacomI (|2006l ) summarise some of these proposals. 
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